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Abstract 

Group theoretical analyses are applied to the magnetic and electronic structures of 
NdNiC>3 and YMO3 , whose electronic structures have been studied very recently by 
the LSDA+U method. The space groups of the crystallographic structure of these 
materials are Pbnm and P1\jn, respectively. [1,2] There are no Jahn- Teller distor- 
tion modes consistent with observed propagation vector of spin polarized neutron- 
diffraction. Two possible spin configurations (magnetic space groups) C a mc2\ and 
P\,a are derived from Pbnm for NdNiO"3. Then the electronic structures with each 
symmetry are studied by using the tight-binding model. Clear separation of partial 
density of states exists in P^a but does not in orthorhombic C a mc2\. The combi- 
nation of results of the present work and the previous LSDA+U calculation shows 
that the ground state magnetic structure of NdNiO"3 is monoclinic P^a. NdNi03 has 
only one kind of Ni sites, while YMO3 has two inequivalent Ni sites. The magnetic 
structure of YMO3 is monoclinic Pf,2i/a, if one assumes no magnetic moment on 
Ni ions (Ni 4+ ) on one of these inequivalent sites. 
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The perovskite nickelates NdNi03 and YNi03 are anti-ferromagnetic insula- 
tors in their low temperature phases. Their magnetic structures are charac- 
terized by a propagation vector k — (|,0, |) [3] and, therefore, the magnetic 
unit cell is a 2 x 1 x 2 supercell of the original crystallographic unit cell. The 
spin alignment in these two nickelates is quite unique one such that the sym- 
metrical coexistence of ferro- and antiferromagnetic couplings along the three 
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pseudo cubic axes as • ■ -Ni T - O -Ni T - O -Ni- 1 — O -Ni- 1 • • -. In YNi0 3 , 
the charge disproportionation 2Ni 3+ — > Ni 2+ +Ni 4+ and the symmetrical co- 
existence of ferro- and antiferromagnetic couplings are both observed, though 
the Ni 4+ ions have small magnetic moments. [2] Diffraction experiments in 
NdNi03 do not show accompanied Jahn- Teller distortion associated with the 
orbital polarizations of Ni e g orbitals. 
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Very recently, we have calculated the electronic structures of the perovskite 
nickelates NdNi0 3 and YNi0 3 by the LSDA+U method which will be pub- 
lished elsewhere. [4] We observed in that study no orbital ordering in NdNi0 3 
and the charge disproportionation in YNi0 3 . The basic consideration by the 
group theory about the magnetic structures, the possibility of the Jahn- Teller 
modes and the simple picture of the e g bands in these structure are given in 
the present report. 

Crystallographic space group: The crystallographic space groups of NdNi0 3 
and YNi0 3 are orthorhombic Pbnm and monoclinic P2 1 /n, respectively. P2\jn 
is the subgroup of Pbnm and their group elements are as follows; Pbnm = 
P2 1 /n®{E, K|±±0}}, P2 1 /n = {E, {a y \\\\}, /, {C 2y \\\\}}. where E is the 
identity operator, / the inversion operator, o~ x the mirror operation with re- 
spect to the y-z plane and C 2y the two-fold rotation operator with respect to 
the y axis. Both NdNi0 3 and YNi0 3 systems contain four Ni ions in their crys- 
tallographic unit cells and each locates on the inversion center. Four Ni sites are 
labeled as Ni(l) = (|, 0, 0), Ni(2) = (0, \, 0), Ni(3) = (|, 0, \) and Ni(4) = (0, \, \). 
In NdNi0 3 , all four nickel positions are equivalent. In YNi0 3 , no symmetry op- 
erations transfer a position Ni(l) or Ni(4) to that of Ni(2) or Ni(3). Therefore, 
there are two inequivalent positions of Ni ions in YNi0 3 . 

Magnetic space group: First we consider the magnetic space groups generated 
from the crystalline space group Pbnm with the experimentally observed prop- 
agation vector k — (|, 0, |). The translation vector along the a-axis of crystal- 
lographic Pbnm structure is denoted by T a . The system is the orthorhombic 
base centered lattice and the (orthogonal) primitive translation vectors should 

be {2T X) T y , 2T Z }. The nonorthogonal translation vectors T x + T z and T x — T z 
transform lattice points at corners of the Bravais lattice to those at the base 
centers. 

Assuming the orthorhombic magnetic unit cell, one of possible spin arrange- 
ments is the structure where Ni(l), Ni(2), Ni(3), Ni(4) in a unit cell of Pbnm 
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have the same spins and the translation operators T x or T z flips the spins of 
Ni ions. [3] Then the ferromagnetically aligned spins are confined within a 
one-dimensional rod extending towards y-direction. The doubled translation 
vectors along the x- and z-axes prohibit the screw operator 2\ along the x- 
and z-axes. [5] Then the unitary part of the resultant magnetic space group 
should be Pmc2 1 = {E, {C 2y \^^}, {cr x \^0}, {a z \00^}}. An augmenting op- 

erator with the time-reversal operator is the translation one T x and the mag- 
netic translation vector should be 2T X . The resultant magnetic space group is 
base centered orthorhombic C a mc2 1 . 

Another possible spin arrangement associated with Pbnm and k = (|,0, |) 
is that where the same spins are located on Ni ions in a parallelepipedon 
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spanned by translation operators {T x , T y , T z — T x \. In this structure, the two- 
fold rotational operator with respect to the z-axis and the mirror operator 
with respect to the x-y plane are prohibited because these contradict to the 
monoclinic Bravais lattice. The screw operator 2± along x-axis is also prohib- 
ited due to cell doubling. There remain only two symmetry operators E and 
the glide operator along the monoclinic axis T z — T x . Because the 
augmenting operator with the time- reversal operator is the translation T x , the 
resultant magnetic space group is monoclinic P^a. The magnetic structure of 
NdNiC>3 with P^a symmetry is of double layers of the (101) plane, on which 
all Ni ions have the same spins. 

In YNiC>3, there are two inequivalent Ni sites, as already stated. If there exists 
no magnetic moment on Ni sites, e.g. Ni(2) and Ni(3), the system has the 
inversion symmetry around the sites of Ni(l) and Ni(4), and the magnetic 
space group should be Pb2i/a. This is actually the case for YNi0 3 where 
the charge disproportionation 2Ni 3+ — > Ni 2+ +Ni 4+ occurs and the magnetic 
moment of Ni 4+ is zero. 

John-Teller distortion: In many transition metal perovskites, the Jahn- Teller 
distortion coexists with the complicated spin and orbital orderings. In per- 
ovskite nickelates, no Jahn- Teller distortion has been reported yet. A pos- 
sibility of the Jahn- Teller mode consistent with the space group Pbnm and 
the propagation vector k = (|,0, |) can be discarded in a following way. A 
unit cell of Pbnm contains four Ni0 6 octahedra which share O ions with each 
other. Let us consider the E g Jahn- Teller modes in two octahedra Ni(l)0 6 
and Ni(4)C>6 which do share no O ions. There are two two-dimensional irre- 
ducible representations in the /c-group G(k = (|,0,|)), and they are listed in 
Ref. [3]. The only elements of non-zero character are E and {o^lHI}} an d 
x({c?,|§f §}}) = ±2. Therefore, the projection operator on those representa- 
tions are P = d/g Y^RaG x{R)*R = §({£|000}}±K|||±}}. When K|±±±}} 
operates on Q 1 u, resultant mode equals to Q A u and vice versa. Here Q^u is 
the E g u mode of Ni(j)06 (j = 1, 4). Then basis vectors of the lattice dis- 
tortion that belong to these irreducible representations are Q l u ± Q 4 u in two 
octahedra Ni(l)C>6 and Ni(4)C>6. The propagation vector k = (|, 0, |) requires 
that the translation operation of T x or T z changes the sign of a displacement 
vector. The resultant displacement vectors are shown in Fig. 1. Now the dis- 
placements of all O ions of Ni(3)Og have been determined, because they are 
contained either in Ni(l)C>6 or in Ni(4)Og. The resultant distortion of the 
Ni(3)0 6 octahedron is an ungerade mode, as seen in Fig. 1, and cannot be 
generated from E g mode on this octahedron. The same arguments can be ap- 
plied to other distortion modes. This is the proof that the E g mode in one 
octahedron cannot be transferred to all octahedra over the whole unit cell. 

Structure of e g bands: Because the 3d e g orbitals in Ni ions are partly filled, the 
highest occupied (HOMO) states should be e g bands. Since the structures of 
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NdNiOs and YNiC-3 have been discussed group theoretically, the band struc- 
ture of the e g orbitals can be analyzed from the similar group theoretical view 
point and this analysis would be useful to understand the calculated elec- 
tronic structure of these systems by the LSDA+U. [4] The discussion can be 
based on a tight-binding Hamiltonian of single spin component and one can 
neglect electron transfer from one spin site to sites with opposite spin. Since 
there are no Jahn- Teller modes and the tilting angle of the NiC>6 octahedra 
is small, distortion can be neglected and an ideal cubic structure is assumed. 
The Hamiltonians do not have any terms of the Slater-Koster parameter V ddlT 
due to this high symmetry. 

First, let us discuss the case of C a mc2 1 . Electrons can transfer only within 
a one-dimensional rod and the ferromagnetic cell extending along the y-axis. 
The e g orbitals are located at four vertices of a unit cube, i.e (0, 0, 0), (0, 0, 1), 
(1,0,0), (1,0,1), with a translation vector t = (1,1,0). The parameter Vdda 
and Vdds are chosen to equal to —1 and —0.1, respectively, and the diagonal 
atomic energy to be zero. Figure 2 shows the resultant energy bands E(k) 
in comparison with those by the LSDA+U method. [4] Because this one- 
dimensional lattice is bipartite, the E{k) curves are symmetric with respect 
to E = 0. As occupation in Ni case is just one electron per each site, the 
Fermi energy is at the energy zero Ep = 0. The result of the tight-binding 
Hamiltonian can reproduce the essential feature of the e g band. Two widely 
dispersive and two narrow bands are occupied and they have no dispersion 
along x- and z-axis. Two dispersive curves are crossing at the Fermi level. The 
k- value of the crossing point actually depends on the ratio Vdds/Vdda, because 
the two-fold degeneracy at the Fermi level is not due to the crystallographic 
symmetry. The separation of partial density of states can be seen very clearly 
in the partial density of states in the present tight binding model, as shown in 
Fig. 2c. The model Hamiltonian keeps much higher symmetry than the actual 
system and the actual system opens the gap as seen in Fig. 2a 

Next, let us consider the case of P^a. The (101) plane of the Pbnm structure 
corresponds to the (111) plane of the cubic lattice. Hereafter, the subscript 
"c" denotes the cubic lattice as (11 l) c . Centers of atomic e g orbitals (Ni ions) 
locate at vertices of a unit cube; (0, 0, 0) and (0, 0, 1) with a translation vector 
a = (1, 0, —1) and b = (0, 1, —1). Since the system has the three- fold rotational 
symmetry around [lll] c axis, the Brillouin zone is a hexagon. The unit cell 
contains two Ni atoms and, therefore, four bands exist in the Brillouin zone. 
Two of them have the dispersion proportional to the Vdds and the other two 
have larger dispersion. They touch each other at the Fermi energy E F . Each 
band is symmetric around Ep, because the lattice is bipartite. The contour plot 
of the band energy E(k) touching Ep is shown in Fig. 3, where the Vdda and 
Vdds are set to be —1.0 and —0.1 respectively, together with the band structures 
of the LSDA+U calculation. [4] The Fermi surface is reduced to just points, six 
vertices of the first Brillouin zone (BZ), and there is two- fold degeneracy at the 
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Fermi surface due to three-fold rotational symmetry. Therefore, this system is 
a marginal metal. The eigen states at the vertex of the BZ are corresponding to 
bonding and anti-bonding states of tp u ± on each adjacent plane, where <p u+ = 
-^(y? 3 z 2 -r 2 + iy?a; 2 -y 2 ) and w_ = ^(</? 3z 2_ r 2 - 'np x 2_ y i). This is not only the 
case at the Fermi energy. The system of the model Hamiltonian is the trigonal 
D,3d and the basis wavefunctions are also tp u ±, which are the basis of the E g 
representation. Therefore, this system does not cause the orbital polarization 
at all, and the partial density of states by the tight-binding model is shown 
in Fig. 3c. Because this system is a marginal metal, the small perturbation 
easily lifts the two-fold degeneracy at the Fermi surface and leads the system 
to the insulating phase. Actually, the gap positions in Fig. 3a are identical to 
the vertices of the BZ in Fig. 3b and the system has no three-fold symmetry. 

Both simple tight-binding models based on C a mc2 1 and P^a spin alignment, 
give metallic bands in contradiction to the fact that NdNiC>3 is insulator, sim- 
ply because the models have much higher symmetry than that of re realistic 
models. The result of LADA+U calculation, which is referred before, shows 
that the bands in P^a open the gap of 0.11 eV without orbital polarization 
and those in C a mc1\ also open the gap of 0.24 eV with a slight orbital or- 
der. [3] Since the orbital polarization has not been observed and no Jahn- Teller 
distortion reported, the spin alignment of P^a is probably the case. 

The charge disproportionation 2Ni 3+ — > Ni 2+ +Ni 4+ is observed experimentally 
in YNiC>3. This is accompanied by the crystallographic symmetry reduction 
from Pbnm to P2\jn. If the spin on Ni 4+ ions is zero, then only one magnetic 
structure Pb2i/a is allowed. In this case, the charge polarization (charge order) 
opens the gap. Our LSDA+U calculation gives the band gap of 1.03 eV in 
YNi0 3 . [4] Experimentally observed magnetic moment on Ni 4+ ions is not 
zero and this discrepancy may be due to the dynamical effect of spins. [2] 

Conclusion: We discussed very details of the magnetic space group for NdNiO"3 
and YNiC>3, and discarded the possibility of the Jahn- Teller distortion extend- 
ing over the whole magnetic cell. The electronic structures of e g bands are 
analyzed by the tight-binding Hamiltonian, and compared with the calculated 
results by the LSDA+U method. [4] 

The present tight-binding model, in both cases of NdNi0 3 and YNi0 3 , is so 
simple and several symmetry elements are not included in the Hamiltonian. 
Then the resultant band structure is of the marginal metal but not of insulator. 
This discrepancy is not serious at all. 

Acknowledgment: This work is supported by Grant-in-Aid for COE Research 
"Spin-Charge-Photon" . 
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Fig. 1. The displacement vectors corresponding to a basis function of group 
of k = (^,0,^) of Pbnm. Center of each octahedron is occupied by Ni ion; 
Ni(l)=(±,0,0), Ni(2)=(0, ±,0), Ni(3)=(±,0, \), Ni(4)=(0, \, \). The Ni ions with 
prime marks can be obtained from the original (without a prime) ions by the prim- 
itive translations. 

Fig. 2. Energy bands E(k) for one spin conponent in NdNiC-3 with C a mc2\ by the 
LSDA+U method (a) and those by the tight binding model (b). The Slater-Koster 
parameters Vddu an d Vdd& eqaual to —1.0 and —0.1 in (b), respectively, (c) The 
partial density of states of e g states in NdNiC-3 by the tight binding model. The 
orbital polarization can be clearly seen. 

Fig. 3. The energy bands E(k) for spin strucrutre P^a of NdNiO"3 by the LSDA+U 
method (a) and those by the tight binding model (b). The results (b) should be 
folded along chain line and the resulting E-k curves correspond to those in (a). The 
Slater-Koster parameter V^dc an d V^ds eqaual to —1.0 and —0.1 in (b), respectively, 
(c) The partial density of states of e g states by the tight binding model. The sepa- 
ration of partial density of states does not exist, (d) The contour plot of E(k) of the 
corresponding tight-binding model Hamiltonian. The g a and gj, are the basis set of 
the reciprocal lattice, and the corresponding translation vectors are a = (1,0, — 1) 
and b = (0, 1, —1). Translation vector c and a corresponding vector g c are parallel 
to (1,1,1). 
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